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We study one-dimensional truncated (no more than 2 particles on a site) bosonic Hubbard model 
in both repulsive and attractive regimes by exact diagonalization and exact worldline Monte Carlo 
simulation. In the commensurate case (one particle per site) we demonstrate that the point of Mott- 
insulator - superfluid transition, (U/t) c = 0.50 ±0.05, is remarkably far from that of the full model. 
In the attractive region we observe the phase transition from one-particle superfluid to two-particle 
ON one. The paring gap demonstrates a linear behavior in the vicinity of the critical point. The critical 

ON ■ state features marginal response to the gauge phase. We argue that the two-particle superfluid is a 

macroscopic analog of a peculiar phase observed earlier in a spin-1 model with axial anisotropy. 

(D ■ 

m ; 

; I. INTRODUCTION 

In the family of Hamiltonians, which the physics of strongly correlated systems deals with, the one-dimensional 
■ bosonic Hubbard one, 

m : 

ON' H = J2 < -t(ala i+ i + H.c.) + -riiim - 1) \ , (1.1) 

i> : 

On . [o| creates a particle on the site i, t is the hopping amplitude, rij = a\a,i, N a is the number af sites] plays an important 
part as a laconic model for the superfluid (SF) - Mott-insulator (MI) phase transitiorap ( and for SF - Bose-glass 
(BG) and BG-MI transitions in the presence of disordeiata). Along with the full model Eq. ( |1.1| ) sometimes actually 
its truncated counterpart, defined by the requirement that site occupation numbers are less than 3, is considered: 
. either deliberately, for the purpose of simpli£gationQ, or implicitly, in approximate treatments which do not distinguish 
between the full and the truncated modelsErtj. One ordinarily motivates such a replacement by the observation that 
in the full model even at the boundary point U — [At U < the system is unstable against a collapse of almost all 
O ■ the particles onto one site.] the quantum-mechanical probability to find out that the occupation number of a given 
site is greater than 2 equals to (e — 5/2)/e « 0.080, so it is very likely that the triply (and more) occupied sites play 
no essential role. _ 

However, a rather accurate method "exact diagonalization plus renormalization group (RG)"tl demonstrated the 
absence of the superfluid phase in the truncated model almost down to U ^ D, in sharp contrast to the above- 
mentioned reasoning, and to the results of a number of approximate treatmentSLTLl. It is essential that the conclusion 
of Reffl about the absence of superfluidity in the macroscopic limit is not a direct one. What was really demonstrated 
is that at arbitrarily small U the known RG flow of superfluid parameters, initialized numerically at intermediate 
scales, leads to the destruction of superfluidity in favor of Mott insulator. Hence, it would be instructive to check this 
fact explicitly by Monte Carlo (MC) simulation of a sufficiently large system, allowing direct measuring the insulating 
gap at U <C 1. This is the primary goal of the present paper. 

Another interesting aspect of the truncated bosonic Hubbard model, distinguishing it from the full one, is that at 
U < it is stable against the collapse onto one site. Therefore our secondary goal is to investigate the region of 
attraction: First, one may expect to find here the superfluid phase. Secondly, at a sufficiently strong attraction a 
pairing of particles should take place, as is clear from the perturbative analysis at U — > — oo, and it is interesting to 
study numerically the corresponding phase transition. Note that in Ref.EI the close problem of the BCS-like pairing 
was examined in a dilute atomic Bose gas with an effective attractive interaction. 

Finally, ID truncated bosonic Hubbard model is not of academic interest only. Macroscopically, it is very close to 
a spin-1 (anisotropic) chain in an easy-axis (U < 0), or easy-plane (U > 0) environment. Its study thus provides deep 
insight into the physics of phase transitions in the spin chains. 

In the present paper we report a study of the model by the exact continuous-time worldline MC approach, developed 
recently by Prokof 'ev, Svistunov and TupitsynE2Hiil. (Some preliminary results were obtained by the standard worldline 
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MC method. For some special purposes we also employ exact diagonalization.) In the commensurate case we have 
found the critical value of the MI-SF transition, (U/t) c = 0.50 ± 0.05, which is more than 6 times lower than that of 
the full model. In the attractive region, U < 0, we observe a specific pairing phase transition when the system passes 
from the one-particle superfluid vacuum to the two-particle one. At the filling factor equal to unity [commensurability 
being not relevant] the transition occurs at U/t w —6.0. The critical behavior is characterized by linear collapse of 
the pairing gap as a function of U , and by a marginal response to the gauge phase at the transition point. 

We also consider spin-1 analog of the truncated bosonic model, which was studied earlier_by exact 
diagonalizatior£j£j|— .The question of our interest here is a peculiar phase which was observed in Ref£3, but was 
not revealed in Ref.O, and, as far as we know, was never discussed since then. We demonstrate that the phase does 
really exist, being a direct macroscopic analog of the two-particle superfluid phase of the bosonic system. 

II. REPULSIVE REGIME 

The MI-SF transition in the truncated one-dimensional commensurate bosonic Hubbard model was studied previ- 
ously by a number of approximate methods. In Ref.uthe real-space renormalization-group method was used, with the 
result (t/U) c = 0.215 for the critical point. In Ref.Q Bethe-ansatz approximation for the full model was considered, 
with a conjecture that the approximation actually corresponds to the exact solution of the truncated model. The 
critical value obtained was (t/U) c — 0.289. In Ref.Q the strong-coupling expansion for the full model was exam- 
ined. The results [(t/U) c = 0.215 for the bare third-order expansion, and (t/U) c — 0.265 with a correction to the 
Kosterlitz-Thouless critical behavior] are quite applicable to the truncated model, since the triply and more occupied 
sites contribute only to the higher-order terms. 

The above-cited results are more or less close to each other, and to the critical value of the full model, (t/U) c = 
0.304 ±0.002, obtained with a controllable accuracy by an accurate method "exact diagonalization + renormalization- 
group (RG) analysis''^. [Practically the same-.value, (t/U) c = 0.30, was found.-also by density-matrix real-space RG 
methodL3, and by accurate MC simulations!!^.] However, the method of Reffl demonstrated that in the truncated 
model the critical value should be essentially shifted towards small {7's: at t/U = 0.3 the truncated system was found 
well inside the insulating region. Moreover, the very fact of the existence of superfluidity at U > was questioned 
[with a reservation that the maximal available size of the cluster, 16 sites, was not enough to resolve unambiguously 
the behavior at small U's.] 

To clarify this situation, we perform an accurate MC study of the commensurate truncated model. We start with 
the standard worldline MC methods and the standard calculation of the insulating gap A = /i + — where the 
chemical potentials fi + and /z_ are given by the relations fi+ = E(N a + 1) — E(N a ) and [i- = E(N a ) — E(N a — 1), 
E(Nb) being the energy of a system consisting of Nb bosons. In \i± vs t/U diagram the region of the insulator phase is 
bounded by the curves fj,+ (t/U) and fi-(t/U). This diagram for N a = Nb = 50 is shown in Fig. 1, in comparison with 
that for the full system. At large U/t phase boundaries in the truncated case are close to those of the full one. This 
is the region, where the strong-coupling approximation is rather accurateu. At smaller U/t, however, the difference 
between the full and truncated models is clearly seen. At U/t < 3.3, where the full model is already in the superfluid 
phase, the truncated model demonstrates a small, but quite pronounced insulating gap. It is seen also, that the 
system is rather subtle in this region, as the gap is considerably smaller than the hopping amplitude. This feature 
makes it difficult to extract the critical point from the data presented, even with the RG-analysis of the gap, since 
the computational errors (~ 0.05i) are comparable with the typical values of the gap (~ O.li) in a rather extensive 
critical region. 

Fortunately, for the problem of pinpointing the critical ratiOp£//t we had an opportunity to take advantage of the 
continuous-time worldline MC approach with "worm" updatetil. The code described in Ref.til directly applies to 
our problem. Working in the grand canonical ensemble ,j— the code allows to extract the macroscopic parameter K 
(main supexfluid characteristic of one-dimensional systemlla) from the histograms for winding numbers and numbers 
of particlesU. Given the values of K for two different system sizes, one easily obtains macroscopic value of K from 
the first integral of thg RG equations, and is able thus to determine whether the system is superfluid or dielectric in 
the macroscopic limitB. The sizes of the systems we used for this purpose were N a — 128 and N a = 16. This way 
we obtained the critical point of the SF-MI transition in the truncated model. (The details of calculation and the 
estimation of the error see in RefH) 

(U/t) c = 0.50 ± 0.05 (2.1) 
As expected, this critical ratio is several times lower than that of the full model. 
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III. ATTRACTIVE REGIME 



Consider the truncated model with the attractive interaction, U < 0. The behavior of the truncated model in this 
region differs qualitatively from that of the full one, the latter being just trivial because of the collapse of almost all 
of the particles onto one site. The constraint on the occupation numbers renders the truncated model stable against 
the collapse, so the point U — is not a special one for it. 

In Fig. 2 we show the \x vs U/t diagram for the truncated model in the wide range of variation of U ft, including 
the attractive region up to | U \ /t = 8. In the region —6.0 < U/t < 0.5 the one-particle gap vanishes, indicating 
that the ground state is an ordinary one-particle superfluid. At U/t < —6.0 the one-particle gap (A ~| Z7 |) appears 
again. Apparently, the point U/t < —6.0 corresponds to the pairing of the carriers. To demonstrate the liquid 
behavior of the two-particle state, we calculate two-particle gap A2 [A2 = /J2+ — ^2-, A*2+ = E(Nb + 2) — E(Nb), 
/i2- = E(Nb) — E(Nb — 2)]. The results are presented in Fig. 3. These reveal the gap corresponding to one-particle 
Mott insulator at U/t > 0.5, and absence of the gap at negative U's. 

In Fig. 4b we show the one-particle gap at the critical region. Within the accuracy of our results, the gap varies 
linearly with U /t. 

Obviously, the pairing transition, being of liquid-liquid type, should take place for non-commensurate fillings as 
well. [We checked this explicitly for certain non-commensurate fillings.] Though not relevant for the pairing transition, 
commensurability may lead however to the formation of the Mott insulator for pairs. The absence of the two-particle 
gap in Fig. 3 at U < means that the state of two-particle Mott insulator does not take place, except probably the 
case of very large | U \ . 

The limit of large, .attraction can be studied analytically. By the standard second-order perturbative treatment in 
the parameter t/[/tM3 one obtains the following effective Hamiltonian in the subspace of states with doubly occupied 
or empty sites: 

t 2 Na 

1 ' i— 1 

where b\ creates a pair on the site i. The states with single occupations of sites are separated from the ground state by 
the gap close to | U |. We thus come to the familiar ID half-filled hard-core model with the parameters corresponding 
just to the critical point of superfluid - insulator (commensurate density wave) transition (see, e.g., El). So we conclude 
that in our model the "critical point" for the Mott transition for pairs is U = 00, and for finite U's the transition 
does not take place. 

IV. RESPONSE TO GAUGE PHASE 

Pairing leads to a dramatic change in the response of the system to the global gauge phase. It is instructive thus 
to study this response in the vicinity of the critical point. 

We introduce the gauge phase, to the model by conventional transformation of the hopping amplitude: t — > 
texp(i<$>/N a ). 

First we study the groundstate energy, Eq, as a function of <&, at different U's, by exact diagonalization for the 
chain N a = iV& = 12, see Fig. 5. [The results for the groundstate persistent current, Jo(&), immediately follow from 
Jo( ( E > ) = dEo(Q)/d&.] Curve 1 corresponds to the well-defined one-particle ground state laying in the zero momentum 
sector at any value of The well-defined two-particle ground state (curve 4) lays in the zero sector while $ varies 
from to 7r/2. At $ = n/2 it jumps to the momentum sector N^/2 (in the units of minimal non-zero momentum), 
and remains there till $ = 37r/2. The curve Eq($>) acquires a characteristic cusp form, symmetric with respect to the 
line <& = 7r/2. Curves 2 and 3, corresponding to the cross-over region, demonstrate a marginal behavior of Eq(&): 
The ground state lays in the zero momentum sector until the phase reaches some critical value (greater than 7r/2), 
and then jumps to the sector Nb/2. The curve is not symmetric with respect to the cusp point. 

We see that the marginal response to the gauge phase in the vicinity of the critical point is due to the competition 
(at <I> > 7r/2) between the groundstate in zero sector and that in the sector Nt>/2. At finite, but sufficiently low 
temperatures both states contribute to the marginal response (provided the system is close enough to the critical 
point), with their Gibbs factors. As an example, we present in Fig. 6 the equilibrium persistent current as a function 
of $ obtained by MC. 

The marginal response to the gauge phase can be used to obtain a rather accurate extrapolation for the critical 
point. At N a — > 00 the size of the cross-over region, where the marginal response takes place, tends to zero. Hence, 
we may consider the value of U, corresponding to some characteristic marginal-response feature, as a function of 
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N a , and extrapolate to N a — oo. For one thing, one may consider just the point of the appearance of the marginal 
response, i.e. the point where at the phase $ = ir the energies in the zero sector and in the sector iV&/2 coincide. 
Corresponding extrapolation from the exact-diagonalization results is presented in Fig. 7. The difference between the 
results of linear and quadratic extrapolations provides a reasonable estimate for the error. 



V. SPIN-1 CHAIN WITH AXIAL ANISOTROPY 



Macroscopically, the truncated bosonic model considered (TBM) in the previous sections is analogous to the spin-1 
chain with axial anisotropy (SI): 

H = -lY,( s ? s 7 + s ^ s l) + v l>2 s * s 3+ u Y,( s tf > 

(ij) (ij) 1 

where U is the parameter of the interaction with the easy-magnetization axis, if U < 0, or easy- magnetization plane, if 



U > 0. [Generally speaking, the term with V renders the model (p.l|) much reacher than the model (1.1). Obviously. 



to make the analogy more close, one should also add a term with nearest-neighbor interaction to the Hamiltonian 



(1.1). To our purposes, however, this is not essential]. . .. . 

The zero-point phase diagram of the model (5.1) was studied in Refs.HLij by exact diagonalization on the chains 



up to N a — 12, see Fig. 8. The phases were identified as follows: 

1) ferromagnetic ((S z ) — ±N a ); 

2) antiferromagnetic ((S z ) — in the ground state as well as in the first excited one); 

3) Haldane gap state (the state with the gap in the spectrum, and (S z ) = ±1 in the first excited state); 

4) XY-sp'm liquid (the gapless state with (5*2) = ±1 in the first excited state). 

[Haldane gap state in SI is equivalent to Mott state in the TBM, while and the XY-sp'm liquid in SI is equivalent 
to one-particle superfluid in the TBM.] 

Besides, in Ref.113 one more phase was observed, the so-called "spin-l/2-like XF-phase" (region 5 in Fig. 8), .-which 
was characterized by (S z ) = in the ground state and (S z ) = ±2 in the first excited state. But later, in Ref.E3 this 
phase was not detected, and, as far as we know, it had not been mentioned in literature since that time. 

Our aim is to make sure that this phase does really exist and is equivalent to the two-particle superfluid in 



TBM. To this end we perform MC simulation of the model (5.1) in the bosonic representation (Holstein-Primakoff 
transformation) for N a = 50, at V = — 0.05i. The fxvsU/t diagram is shown in Fig. 9. It demonstrates the remarkable 
similarity with the behavior of TBM near the point of the pairing phase transition: the one-particle gap, which shrinks 
linearly in the wide interval of the parameter U. The critical point (U /t) c rj —2.0 coincides with the value found in 
RefH 



Note that in the limit U — > 00 the model (5.1) with V = also approaches the effective Hamiltonian (p.l|), and the 
equivalence to TBM becomes exact. 
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FIGURE CAPTIONS 

Fig. 1. fj, vs U/t diagram for the truncated model (o) with N a = Nb = 50 obtained-by the Monte Carlo method. 
T = 0.0625t. For comparison the phase diagram of the full boson model (+) is shownlia. 

Fig. 2. /i vs U/t diagram, obtained for N a = Nb = 50 by Monte Carlo in the wide parameter region. 

Fig. 3. fi2 vs U/t diagram for N a — Nb — 14 (exact diagonalization) . 

Fig. 4. One-particle gap at the critical region T = t/50, N a = 50 (Monte Carlo). 

Fig. 5. Groundstate energy vs gauge phase at N a = Nb = 12 (exact diagonalization): (1) U/t = —5.5, (2) U/t = —6.0, 
(3) U/t= -6.5, (4) U/t = -10.0. 

Fig. 6. Equilibrium persistent current vs gauge phase in the marginal region. U/t — —5.92, N a = 50, T = t/150. 

Fig. 7. Scaling for the point of the appearance of the marginal response to the gauge phase for chains with iV a = 
8, 10, 12, 14, 16 (exact diagonalization). 



Fig. 8. Phase diagram for spin-1 chain with axial anisotropy, from Ref.olij. (1) ferromagnetic state, (2) antiferro- 
magnetic state, (3) Haldane gap state, (4) XY-sp'm liquid, (5) "spin-l/2-like ATY-phase". 

Fig. 9. n vs U/t diagram for the spin chain. V/t = -0.05, N a — 50, T — t/50. 
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